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Abstract 


A  description  for  the  strain-rate  and  temperature-dependent  behavior  of  pure  tantalum  (Ta) 
at  large  strains  is  developed.  An  integral  part  of  the  model  incorporates  the  kinematics  of 
crystallographic  slip,  and  thus,  the  rotation  of  single  crystals  within  the  material,  so  as  to  reflect 
the  evolution  of  anisotropy  resulting  from  applied  mechanical  deformation.  Such  deformation 
is  accommodated  via  bulk  dislocation  motion  and  governed  by  interactions  that  may  or  may  not 
be  thermally  assisted.  The  model  represents  each  discrete  slip  system  as  a  single  facet  in  a 
multisurface  plasticity  theory,  which  is  well  suited  to  high-rate  numerical  methods  (explicit 
integration  schemes).  A  formulation  of  this  type  allows  for  the  complete  kinematic 
decomposition  of  macroscopic  material  rotations  and  the  rotations  of  single  crystals  due  to 
motion  through  the  lattice.  Applications  to  high-rate  deformation  in  polycrystals  is  conducted 
within  a  finite  element  implementation  by  invoking  a  Taylor  criterion  for  interaction  between 
the  macroscopic  and  the  mesoscopic  (single-crystal)  length  scales.  The  model  behavior  is 
examined  in  application  to  high-rate  problems  with  increasingly  complex  geometries 
(homogeneous  uniaxial  compression,  the  impact  of  a  textured  cylindrical  bar  into  a  rigid  anvil, 
and  the  explosive  deep-drawing  of  a  Ta  disk).  The  final  example  compares  through-thickness 
deformation  in  the  Ta  disk  as  a  Action  of  initial  texture  and,  therefore,  emphasizes  the  profound 
effects  of  plastic  anisotropy  during  finite  deformation. 
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1.  Introduction 


By  its  very  nature,  crystallographic  glide  on  a  restricted  set  of  slip  directions  produces 
anisotropic  plastic  behavior  in  the  single-crystal  constituents  of  a  polycrystalline  aggregate. 
Even  though  a  random-oriented  polycrystal  may  exhibit  isotropic  mechanical  properties,  fi¬ 
nite  deformation  in  polycr3rstalline  aggregates  is  often  accompanied  by  severe  reorganization 
of  the  material’s  microstructure.  Upon  such  reorganization,  initially  random  distributions 
of  single-crystal  constituents  may  begin  to  align  in  preferred  orientations  (texture).  There¬ 
fore,  within  materials  that  deform  by  crystallographic  slip,  dislocation  motion  causes  grain 
rotations  and,  thus,  after  large  deformation,  may  cause  severe  mechanical  anisotropy  even 
in  materials  that  are  initially  isotropic. 

Such  large  strains  may  arise  within  certain  general  engineering  applications,  but  more 
often  they  may  arise  during  the  thermo-mechanical  processing  of  stock  material  used  in  the 
production  of  engineering  systems.  In  the  latter  case,  analysis  of  these  systems  may  be  crit¬ 
ically  dependent  on  knowledge  of  the  initial  material  state  and  inherent  material  anisotropy. 
For  example,  the  construction  of  vehicle  body  panels  from  rolled  sheet  involves  a  material 
with  an  initially  high  degree  of  anisotropy  due  to  the  previous  rolling  procedure.  An  analysis 
of  such  a  structure  may  have  to  account  for  initial  anisotropy,  but  the  evolution  of  anisotropy 
during  the  life  of  the  structure  would  not  be  great.  High-strain-rate  applications,  however, 
often  involve  large  deformation  of  material  that  has  already  undergone  extreme  microstruc- 
tural  reorientations  during  the  initial  processing  (i.e.,  manufacturing  process)  and  now  faces 
further  large  deformation,  microstructural  evolution,  and,  thus,  evolving  anisotropy  during 
its  performance  lifetime.  More  importantly,  evolving  anisotropy  may  be  critically  important 
in  setting  the  stage  for  late-strain  localization  and  ultimate  structural  failure.  Accurate  anal¬ 
ysis  of  such  devices  depends  on  understanding  how  material  anisotropy  evolves  during  the 
deformation  of  complex  geometries  and  how  such  anisotropies  affect  the  dynamic  evolution 
of  macroscopic  features. 

Macroscopic  mechanical  anisotropy  evolves  with  increasing  deformation  via  a  number  of 
mechanisms.  In  brittle  materials,  changes  in  the  preferred  deformation  behavior  is  usually 
a  function  of  damage  on  either  the  micro  or  macro  scale.  In  ductile  metals,  such  as  Ta, 
directional  dependence  evolves  via  several  mechanisms  that  commonly  involve  the  motion 
of  dislocations  through  the  crystallographic  lattice  of  single  crystals,  and  in  polycrystalline 
materials,  the  interaction  and  reorientation  of  these  crystals. 
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Observations  of  macroscopic  phenomenology  during  the  deformation  of  anisotropic  mate¬ 
rials  have  led  to  the  formulation  of  anisotropic  yield  surfaces  [1,  2]  and  kinematic  hardening 
models  [3]  that  are  capable  of  reflecting  the  macroscopic  directionality  of  such  materials. 
These  models,  however,  are  only  a  description  of  observations  and,  even  when  well  cali¬ 
brated,  will  never  become  predictive  at  large  deformation  when  the  initial  anisotropy  begins 
to  evolve.  Within  this  framework.  Maudlin  et  al.  [4,  5]  have  sought  to  incorporate  the 
crystallographic  orientations  into  an  anisotropy  model  by  direct  measurement  using  x-ray 
diffraction.  These  measured  orientations  are  then  used  to  formulate  a  faceted  yield  surface  as 
a  convex  flow  envelope  representing  the  initial  grain  orientations  in  the  material.  Although 
ideally  suited  for  high-rate  applications  and  physically  based,  this  yield  description  lacks 
evolution  of  grain  orientations  at  large  strains. 

Physically  based  models  for  texture  evolution  have  been  under  development  for  a  number 
of  years  [6,  7]  but  have  only  recently  been  applied  to  deformation  at  high  strain  rates  [8,  9]. 
Applications  of  these  models  have  mainly  been  concerned  with  the  analysis  of  homogeneous 
deformations.  Robust  polycrystal  implementations  for  the  analysis  of  nonhomogeneous  de¬ 
formation  fields  are  indeed  quite  new  [10,  11].  Most  recently,  Kothari  et  al.  [12]  have  de¬ 
scribed  an  elasto-viscoplastic  constitutive  model  for  polycrystalline  Ta  deforming  at  high 
strain  rates.  This  model  has  been  implemented  into  the  finite  element  code  ABAQUS  [13] 
and  used  to  analyze  the  impact  of  a  pretextured  Ta  rod. 

Appropriately,  this  work  is  concerned  with  describing  how  Ta  and  other  ductile  poly¬ 
crystals  behave  when  deforming  to  large  strains  at  high  rates.  We  begin  with  a  description 
of  the  fundamental  kinematics  of  crystallographic  slip  in  single  crystals,  incorporating  a 
likely  constitutive  description  for  the  elasto-viscoplastic  behavior  of  single  crystals  and  an 
appropriate  numerical  implementation  of  the  combined  constitutive  theory.  This  single¬ 
crystal  behavior  is  used  in  an  appropriate  averaging  technique  to  achieve  a  description  for 
polycrystalline  material  behavior  that  reflects  the  evolution  of  material  anisotropy  via  the 
underlying  physical  mechanisms.  This  model  is  especially  suited  for  implementation  into 
a  coupled  thermomechanical,  large-strain,  high-rate  finite  element  formulation  for  solving 
complex  boundary  value  problems.  Example  calculations  include  both  homogeneous  and 
nonhomogeneous  uniaxial  compression  of  random  and  pretextured  Ta  stock.  The  model  is 
then  used  to  investigate  through-thickness  sheet  thinning  during  explosive  deep-drawing,  a 
dramatic  example  of  texture-induced  plastic-anisotropy  with  practical  implications. 
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Notation  used  in  this  work  is  based  on  the  following  convention.  Scalars  are  denoted 
by  mathematical  italics,  {A,  a,  a),  vectors  by  lowercase  boldface  (a),  second-order  tensors 
by  uppercase  boldface  (A),  and  fourth-order  tensors  by  calligraphic  uppercase  boldface  (K). 
The  summation  convention  is  used  for  Latin  indices.  Superscripts  pertaining  to  slip  systems 
are  written  in  Greek  letters;  summation  over  these  are  indicated  explicitly.  Dot  and  double 
dots  are  used  to  indicate  the  following  products: 

F  •  m  ^  Fijirij,  C  :  C  ^ 


and 

/C  :  L  ^  ICijkiLik',  T  •  D  ^  TikDkj  • 
Finally,  time  derivatives  are  denoted  by  superposed  dots,  e.g.: 

F  =  dF/dt, 

while  the  transpose  of  a  tensor  is  designated  by  the  superscript  T. 


2.  Material  Description  for  Crystallographic  Slip 

The  following  section,  beginning  with  a  description  for  the  thermally  activated  glide  of 
dislocations  along  a  restricted  set  of  slip  directions  and  plane  normals  followed  by  the  kine¬ 
matic  implications  of  restricted  glide,  constitutes  the  current  model  for  thermo-mechanical 
deformation  in  a  single  crystal. 

2.1  Description  for  Thermally  Activated  Dislocation  Glide 

A  framework  for  the  bulk  motion  of  dislocation  glide  that  is  both  temperature  and  strain 
rate  dependent  [14, 15]  has  been  successfully  calibrated  to  Ta  and  Ta-W  alloys  under  a  broad 
range  of  temperatures  and  strain  rates  [16].  Most  recently,  Kothari  and  Anand  [12]  incor¬ 
porated  a  similar  model  to  explicitly  describe  crystallographic  slip  on  discrete  slip  systems 
and  slip  planes.  The  following  section  is  a  brief  review  of  how  dislocations  may  interact  with 
various  obstacles,  and,  how  these  obstacles  may  be  affected  by  changes  in  temperature  and 
strain  rate. 

Dislocation  motion  is  assumed  to  result  from  both  applied  stress  and  statistical  fluc¬ 
tuations  due  to  thermal  vibrations.  The  motion  is  affected  by  obstacles,  both  long  and 
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short  range  with  respect  to  lattice  dimensions;  long  range  obstacles  typically  include  grain 
boundaries  or  large,  second-phase  particles  while  the  short  range  barriers  may  include  forest 
dislocations,  interstitial  or  solute  atoms,  small  precipitates,  or  in  the  case  of  pure  Ta,  large 
Peierls  stresses  [16,  17,  18].  Such  obstacles  are  said  to  characterize  the  material  “state,”  or, 
in  terms  of  stress,  where  r“  is  the  resistance  to  dislocation  motion  on  the  a  system,  we  say 

r°=T;  +  '£Si{0,\r\)rr,  w 

i=l 

where  r“  is  the  resistance  associated  with  long-range  interactions  that  cannot  be  thermally 
assisted.  Each  type  of  thermally  activated  obstacle  (summed  over  i)  that  may  resist  disloca¬ 
tion  motion  is  characterized  by  r“,  the  stress  required  to  overcome  the  barrier,  or  constant 
“structure,”  at  0  K.  The  temperature  and  strain-rate  contributions  to  overcoming  such  bar¬ 
riers  are  reflected  in  the  Si{6,  |7®|)  terms,  where  9  denotes  temperature  and  7“  denotes  the 
slip  rate  on  the  a’th  slip  plane.  Since  rf  defines  the  maximum  barrier  for  a  given  structure, 
the  Si  are  by  definition  <1.  In  the  previous  equation,  thermally  activated  obstacles  are  con¬ 
sidered  on  the  scale  of  the  crystallographic  lattice  dimensions;  obstacles  such  as  the  Peierls 
resistance  to  glide  may  be  overcome  by  thermal  activation  energy  alone.  Thus,  the  functions 
Si(9,  |7“|)  lower  the  resistance  due  to  the  material  state  Tf  by  accounting  for  the  statistical 
frequency  that  such  obstacles  may  be  overcome  by  thermal  fluctuations.  Each  Si{9,  |7“|) 
describes  the  short-range  interaction  kinetics  via  the  phenomenological  [15]  form 


5(0,iri)i=  1- 


k9\ni^^)Y 


where  b  is  the  lattice  Burgers  vector,  g'o.i  is  the  normalized  activation  energy  for  disloca¬ 
tion/obstacle  interaction,  pi  and  Qi  are  factors  that  represent  geometry  of  i^*  obstacle  type, 
and  7o,i  is  a  reference  shear  rate  taken  to  be  the  same  on  all  slip  systems  [12].  The  dislocation 
kinetics  of  such  flow  is  discussed  in  detail  in  other  works  [14,  15,  19]. 


Kothari  and  Anand  [12],  when  considering  pure  Ta,  allow  only  a  single  thermally  assisted 
obstacle  (representing  the  Peierls  barrier),  and  allow  only  the  athermal  structure  to  evolve 
with  dislocation  motion. 


In  order  to  model  deformation  behavior  in  pure  Ta  over  a  broad  range  of  temperatures 
and  strain  rates,  Chen  and  Gray  [16]  distinguish  between  the  thermally  activated  Peierls 
barrier  and  structure  associated  with  dislocation  hardening  and  recovery  (dislocation  accu¬ 
mulation  and  annihilation).  Further,  they  diflFer  from  Kothari  and  Anand  [12]  by  assuming 
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that  this  second  mechanism  is  thermally  activated.  Following  Chen  and  Gray  [16],  equa¬ 
tion  (1)  becomes 

-1-  SyT^,  (3) 

where  Sp  and  r®  are  associated  with  the  Peierls  stresses  and  do  not  evolve  during  deforma¬ 
tion.  The  Sy  and  r®  pair  are  associated  with  dislocation  density.  Sy  will  take  the  form  of 
equation  (2),  and  r“  will  evolve  according  to  the  Peirce  et  al.  [20]  structure  that  has  been 
applied  within  this  framework  by  Kothari  and  Anand  [12],  i.e.. 


p 

The  hardening  moduli,  h®^,  is  given  as 


(4) 


=[q+{l-q)S^^]h^,  (5) 

where  denotes  the  self-hardening  rate  and  q  denotes  the  ratio  of  self  versus  latent  harding. 
As  in  Chen  and  Gray  [16],  the  self-hardening  due  to  dislocation  motion  will  be  assigned  the 
form  proposed  by  Follansbee  and  Kocks  [15] 


h^  =  h^ 


tanh 

k^] 

S  7-“ 

tan 

MCI 

(6) 


where  r®  takes  the  rate-  and  temperature-dependent  form 


(7) 


In  the  previous  two  equations,  ^  and  r®,o  axe  material  constants,  b  is  the  magnitude  of  the 
Burgers  vector,  k  is  the  Boltzmann  constant,  and  is  a  reference  rate  of  shearing. 


2.2  Single- Crystal  Kinematics 

The  motion  of  an  elasto-plastic  single  crystal  can  be  partitioned  as  the  result  of  elastic 
stretching  accompanied  by  rigid-body  rotation  and  simple  shearing  along  the  various  slip 
systems  of  the  crystal.  The  shearing  motion  is  due  to  dislocation  movement  along  the  slip 
systems  but  is  described  in  a  bulk  continuum  framework.  The  motion  of  the  single  crystal! 
then  can  be  decomposed  into  a  plastic  part,  due  to  crystallographic  glide,  and  an  elastic 
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part,  due  to  lattice  stretching  and  rotation.  Using  F  to  designate  the  deformation  gradient 
of  the  single  crystal,  a  multiplicative  decomposition  [21]  is  written  in  component  form  as 

F  =  F*-FP,  (8) 

where  F^  is  due  to  the  net  motion  of  dislocations  and  F*  contains  the  elastic  distortion  and 
rigid-body  rotation  of  the  lattice  [22,  23,  24].  The  slip  direction  of  a  particular  system,  a, 
is  denoted  by  and  must  lie  in  the  slip  plane  with  normal  The  orthonormal  pair 

s^®)  and  m^®^  convect  with  the  lattice  to  become 

s*(®)  =  F*-  s^®)  (9) 

and 

m‘(®)  =  m(®)-  F*-^  .  (10) 

The  total  velocity  gradient,  L,  can  then  be  written  in  an  additive  decomposition,  i.e., 

L  =  F  •  F-^  =  L*  -h  U,  (11) 

with  L*  and  representing  the  elastic  and  plastic  parts,  respectively.  The  inelastic  defor¬ 
mation  due  crystallographic  slip  is  the  linear  summation  of  slip  along  some  active  subset  of 
the  total  number  of  restricted  glide  systems.  If  we  denote  Mact  l-o  be  this  subset,  then 

Mact  ^act 

L"  =  £  s*(“>  (8>  m*(®)7(®)  =  2  .  (12) 

a  a 

The  plastic  part  of  the  velocity  gradient  can  be  further  decomposed  into  the  symmetric  rate 
of  stretching,  D^,  and  the  antisymmetric  plastic  spin,  QP, 


LP  =  DP  -1-  . 

(13) 

We  denote  the  symmetric  and  antisymmetric  parts  of  the  Schmid  tensor,  P®, 
W®,  respectively,  and  write 

by  R“  and 

r(“)  =  i[s*(®)m*(®)  +  m*(®^s*(®)] 

2'-  ^ 

(14) 

and 

1 

=  i[s*(“^m*(®)  -  m*(®)s*^®^] 

2^ 

(15) 

so  that 

Mact 

DP  =  2 

(16) 

a 
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and 


(17) 


M^act 

=z  . 

a 

Du3dng  deformation,  any  crystallographic  axis,  (e.g.,  the  c  axis)  will  rotate  with  the 
lattice  spin,  Q  as 

c  =  Cl*-  c,  (18) 

where  Cl*  is  the  diflFerence  between  the  total  spin,  fi,  and  the  previously  defined  plastic  spin, 
i.e.. 

Cl*  Cl  -  CV  .  (19) 

3.  Numerical  Implementation 

Finite  element  formulations  for  the  analysis  of  wave  propagation  and  dynamic  inertial  ef¬ 
fects  typically  implement  the  equations  of  motion  in  an  explicit  formulation  where  time  steps 
are  limited  to  some  fraction  of  the  Courant  stability  limit.  As  such,  it  is  common  practice  to 
assume  constant  strain  rates  (as  calculated  from  the  previous  time  step)  across  the  time  step 
and  implement  rate-dependent  material  behavior  as  step-wise  rate  independent  [25].  Stress 
rotations  are  updated  at  half-time  increments  on  either  side  of  the  constitutive  update,  thus 
creating  an  objective  stress  update  decoupled  from  the  constitutive  calculation. 

3.1  Single-Crystal  Return  Algorithm 

Denoting  <r  as  the  rotated,  deviatoric  Cauchy  stress,  the  resolved  shear  stress  on  any 
particular  glide  system,  denoted  by  a,  is  given  as 

t“  =  P“  :  a  .  (20) 

Assuming  that  slip  occurs  when  the  resolved  shear  stress  reaches  some  threshold  value, 
=  r“,  we  can  use  the  previous  equation  to  define  a  yield  function,  for  the  slip 
system  as 

r  =  |*r:P“|-r“  =  0,  (21) 

where  the  components  of  P"  are  the  previously  defined  direction  cosines  that  resolve  <r  into 
the  direction  of  slip.  Note  here  that  for  the  purposes  of  a  numerical  algorithm,  if  r“  <  —if 
(i.e.,  the  resolved  shear  stress  exceeds  the  critically  resolved  shear  stress  so  as  to  create  slip 
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in  the  direction),  then  the  slip  vector,  s*(“\  is  temporarily  changed  to  so  that 

only  positive  7’s  will  be  considered. 

The  collection  of  linear  yield  functions  described  by  equation  (21)  resembles  a  Tresca- 
type  3deld  condition  that  was  actually  described  for  single  crystals  by  Koiter  [26].  Simo, 
Kennedy,  and  Govindjee  [27]  described  a  full  plasticity  theory  for  a  collection  of  smooth, 
multifaceted  yield  surfaces.  The  following  is  an  application  of  the  rate-independent  Simo, 
Kennedy,  and  Govindjee  [27]  theory  applied  to  the  Koiter  [26]  yield  conditions  and  the 
single-crystal  kinematics  described  previously. 

Begin  with  a  trial  elastic  update  of  the  stress  using  a  generalized  form  of  Hooke’s  Law, 

^  ^dd  (22) 


where  C  is  the  material  elastic  moduli  and  At  is  the  integration  time  step  size.  The  calculated 
a  is  then  used  to  evaluate  /“  for  each  slip  system.  If  no  facets  are  violated  (i.e.,  /“  <  0  for 
all  a)  then  <t*'‘  is  accepted  as  the  current  stress  state,  7“  =  0  for  all  a,  and  the  calculation 
continues.  If  the  trial  stress  state  causes  any  /®  >  0,  then  that  facet  is  added  to  a  list  of 
potentially  active  slip  systems.  The  discussion  that  follows  will  be  applied  to  the  set  of  Mact 
potentially  active  yield  facets. 


The  total  plastic  flow  results  from  the  linear  sum  of  plastic  flow  along  each  planar  facet. 
Applying  the  flow  rule  in  this  manner  over  the  Mact  smooth  facets  recovers  equation  (12)  as 


Mact  Mact 

LP  =  V  =  E  . 

^=1  p=i 


(23) 


It  is  important  to  note  here  that  the  scalar  multipliers,  A^,  in  the  plasticity  theory  become 
exiictly  the  slip  rates,  7^,  from  the  single-crystal  theory.  This  allows  for  the  complete  solution 
of  the  coupled  single-crystal/multifaceted  flow  theory.  ' 


Assuming  that  r“  remains  constant  over  the  time  step,  then  the  stress  rate  is  given  by 


=  £  :  L*  =  £  :  {L  -  LP} 


(24) 


and  the  condition  for  plastic  consistency  is  met  by 


(25) 
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Substitution  of  (21),  (24),  and  (23)  into  the  previous  equation  gives  an  Mact  by  Mact 
system  of  equations  for  the  slip  rates  as 

Mact  , 

P“  :  £  :  L  -  2  P“  :  £  :  P^A^  =  0  (26) 

for  a  =  1  to  Mact-  As  discussed  previously,  the  set  of  Mact  systems  is  considered  to  be 
only  potentially  active.  It  is  actually  the  condition  7^  >  0  (not  >  0)  that  activates  the 
slip  system.  Therefore,  each  calculated  7^  must  be  checked  for  activation.  If  7^  <  0, 
then  that  system  is  dropped  from  the  list  of  potentially  active  systems  and  equation(26)  is 
reconstructed  for  the  Mad-1  slip  systems.  The  process  is  repeated  until  all  7^  >  0. 

Once  the  slip  rates  are  known,  the  stress  state  can  be  relaxed  back  onto  the  3deld  surface 

Mo. 

^  -  53  £  :  P^7^At,  (27) 

p=i 

and  the  grain  rotation  can  be  updated  consistent  with  the  single-crystal  kinematics.  If  we 
denote  #  as  the  the  lattice  orientation  of  a  particular  grain  with  respect  to  the  lab  axis,  then 
the  lattice  spin  is  determined  by  Q*  as  per  Asaro  and  Rice  [28]. 

#(t  + At)  =#•#(<),  (28) 


where 

#  =  exp(n*At).  (29) 

We  employ  the  Cayley-Hamilton  theorem  in  order  to  expand  exp(f2)  with  Cl  —  n*At  as 

(30) 


^  sin(n;)-  l-cos(tt;)4.2 
exp(f2)  =1-1-  — —n  + - 


w 


W*‘ 


and  =  1^23  +  ^13  +  ^12- 

As  an  illustrative  example,  consider  the  {110}  <  111  >  family  of  slip  systems  (later,  both 
{110}  <  111  >  and  {112}  <  111  >  slip  are  considered  in  order  to  model  Ta)  under  planar 
flow  conditions.  The  yield  surface  will  be  composed  of  the  inner  envelope  of  hyperplanes 
resulting  from  the  yield  facets  of  each  individual  slip  system.  For  an  incompressible  material, 
the  plane-strain  yield  surface  will  resemble  a  closed  polygon,  which  can  be  expressed  as  a 
single  curve  on  the  {an  -  022)  c^i2  stress  axes  [30].  Such  a  configuration  is  shown  in 
Figure  1. 
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Simo,  Kennedy  and  Govindjee  [27]  give  a  geometric  interpretation  of  this  plasticity 
theory  as  providing  a  closest-point  return  algorithm  to  the  yield  surface.  This  being  the 
case,  we  can  observe  that  the  plastic  straining  will  be  coaxial  with  the  trial  elastic  stress 
only  when  the  trial  stress  is  perpendicular  to  a  particular  yield  facet.  For  example,  when 
aYi  —  —<^22  >  stress  state  will  be  on  the  (110)  [111]  facet  at  the  (an  —  a-i^j^ 

axis.  When  this  is  not  the  case  (for  example,  when  ,  as  shown  in  Figure  1) 

the  resulting  stress  state,  <ti,  will  not  be  coaxial  with  the  straining  increment,  and  after 
successive  straining  increments,  the  resulting  stress  state  will  migrate  toward  a  vertex  in  the 
yield  surface.  Once  the  trial  state  falls  within  a  zone  in  stress  space  bounded  by  normals 
to  more  than  one  yield  facet  (for  example,  a^)  multiple  yield  facets  will  be  active,  and  the 
resulting  stress  state  will  remain  at  the  vertex,  <r2,  until  a  change  in  strain  path  is  significant 
enough  to  move  the  response  away  from  the  vertex. 

Since  the  current  hyperplane  representation  of  each  facet  does  not  actually  incorporate 
any  explicit  description  of  the  yield  vertex,  it  is  possible  for  the  algorithm  to  occasionally 
“mistake”  some  single-facet  stress  states  (trial  stress  states  that  should  only  excite  one  yield 
facet)  for  vertex  stress  states.  Figure  2  shows  two  hyperplanes  from  Figure  1;  it  is,  essentially, 
an  exploded  view  of  Figure  1.  The  trial  stress  state  depicted  in  this  figure  can  easily  be  seen 
to  violate  both  yield  facets  shown.  Of  all  the  possible  stress  states  along  the  yield  surface,  it 
is  actually  <t3,  and  not  the  vertex,  that  minimizes  plastic  work  (closest  point  projection).  It 
is  necessary,  therefore,  to  check  the  vertex  solution  given  by  this  algorithm  against  a  normal 
return  to  each  yield  facet  and  accept  the  solution  that  minimizes  the  plastic  work.  At  this 
point  in  the  algorithm,  the  solution  is  also  checked  for  numerical  accuracy  in  order  to  ensure 
that 

Mact 

E  l/“l  :£  Pl) 

Cl 

where  tol  =  Tc  x  10“^.  If  this  condition  is  violated,  then  the  single  facet  solution  that 
minimizes  the  plastic  work  is  accepted  as  the  final  stress  state. 

3.2  Length- Scale  Interaction 

The  modified  Taylor  polycrystal  averaging  scheme  (Asaro  and  Needleman  [6])  is  incor¬ 
porated  in  order  to  (1)  interact  the  macroscopic  deformation  with  the  local  grains,  and 
(2)  average  the  resulting  stresses  from  the  local  grains  to  the  associated  aggregate  stresses. 
The  model  assumes  that  deformation  within  a  grain  is  uniform  and  that  each  grain  imder- 
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Figure  1.  Six  facets  composing  the  2D  plane  strain  yield  surface  for  the  {110} 
<  111  >  slip  family  in  a  single  crystal  aligned  with  the  macroscopic 
axes  of  stress 
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Figure  2.  Enlarged  view  of  two  yield  facets  from  Figure  1  showing  an  anomaly 
in  the  algorithm  and  the  correct  returned  position 
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goes  the  same  deformation  as  the  entire  aggregate.  Therefore,  the  local  velocity  gradient,  L, 
is  identical  to  the  macroscopic  velocity  gradient,  L.  In  this  way,  deformation  compatibility 
between  the  grains  and  macroscopic  deformation  is  maintained.  After  rotating  the  applied 
L  to  L,  i.e., 

1,9  -  .  L  •  (32) 

the  local  Cauchy  stresses,  a,  are  then  calculated  and  volume  averaged  to  achieve  the  macro¬ 
scopic  stress,  or.  Since  the  calculations  assume  that  all  the  grains  occupy  an  equal  volume 
within  the  polycrystal, 

a  =  —f]  •  o-(»)  •  (33) 

where  K  is  the  number  of  grains  in  the  aggregate. 


4.  Application  to  Dynamic  Deformation  in  Pure  Ta 

The  purpose  of  this  section  is  to  show  how  the  previously  described  model  can  be  in¬ 
corporated  to  resolve  problems  in  dynamic,  nonhomogeneous  texture  evolution;  Increasingly 
complex  examples  of  high-rate  texture  evolution  in  pure  Ta  are  discussed.  In  this  section, 
Ta,  a  body-centered-cubic  (bcc)  lattice  structure  material,  was  taken  to  deform  by  restricted 
glide  on  a  total  of  24  {110}<  111  >  and  {112}<  111  >  slip  systems.  Both  the  slip  resistance 
and  the  evolution  of  slip  resistance  on  each  system  are  taken  as  equal  (i.e.,  Taylor  hardening), 
and  the  material  constants  used  in  describing  the  resistance  to  slip  on  each  slip  system  were 
obtained  from  the  data  reported  in  Chen  and  Gray  [16]  and  scaled  by  a  polycrystal  Taylor 
factor  of  3.58,  as  reported  by  Maudlin  et  al.  [5]  for  the  same  Ta  plate  stock. 

4.1  Uniaxial  Compression 

The  first  example  is  to  model  homogeneous  uniaxial  compression  of  a  Ta  cylinder.  In 
this  example,  a  single  finite-element  was  compressed  under  uniaxial  strain  at  a  constant  rate 
of  10^  to  a  final  true  strain  of  1.0.  The  element  initially  consisted  of  300  random-oriented 
grains;  after  a  true  strain  of  1.0,  the  resulting  texture  is  shown  in  Figure  3.  This  result  is 
consistent  with  high-rate  texture  data  obtained  by  Lee  et  al.  [9]  and  generally  consistent 
with  textures  obtained  in  bcc  materials  (Schoenfeld  et  al.  [8];  Kad  and  Schoenfeld  [29]). 


13 


{Ill) 


(100)  (110) 

Figure  3.  {100},  {110},  and  {111}  pole  figures  for  Ta  compressed  uniaxially  to  a 
true  strain  of  1.0 

4.2  Taylor  Impact  Cylinder 

Maudlin  et  al.  [4]  removed  Ta  cylinder  specimens  from  “heavily”  rolled  Ta  plate.  The 
length-to-diameter  ratio  (L/D)  of  the  specimens  was  5,  and  they  were  removed  such  that 
the  cylinder  axis  was  parallel  to  the  confined  (i.e.,  transverse)  direction  of  the  rolled  plate. 
The  authors  then  launched  the  specimens  from  a  .30-cal.  Mann  powder  gun  at  a  velocity  of 
175  m/s  into  target  4340  steel. 

In  order  to  model  the  previously  mentioned  experiment,  a  finite  element  model  similar  to 
that  used  by  the  authors  was  constructed  using  the  EPIC96  computer  code  and  implemented 
with  the  current  model.  A  total  of  12,960  tetrahedral  elements  were  arranged  in  an  L/D  = 
5  cylinder  and  impacted  into  a  frictionless  surface.  It  is  known  that  the  original  plate  was 
“heavily”  rolled,  so  an  assumed  initial  texture  corresponding  to  plane  strain  compression 
was  imparted  to  48  grains  at  the  gauss  point  of  each  element.  This  initial  texture  is  shown 
at  time  =  0  in  Figure  4,  which  contains  {100},  {110},  and  {111}  pole  figures  from  the  center 
of  the  impacted  surface  before,  during,  and  after  the  forming  event.  From  the  figure,  texture 
evolution  from  the  initially  rolled  material  to  that  characteristic  of  uniaxial  compression 
(Figure  3)  can  be  observed.  Note  that  by  5  fis,  the  initial  texture  has  almost  completely 
disappeared  and  that  the  strong  components  of  the  uniaxial  texture  do  not  begin  to  form 
until  about  50  fjs.  In  spite  of  this  rapid  evolution  toward  an  axisymmetric  texture  at  the 
impacted  surface,  material  away  from  the  impacted  surface  undergoes  a  much  less  drastic 
evolution  and  thus,  retains  a  great  deal  of  the  initial  anisotropy. 

Figure  5(a)  shows  contours  of  effective  plastic  strain  in  the  ovalized  footprint  on  the 
impact  surface.  Even  though  the  texture  at  the  Center  of  the  footprint  evolves  in  an  axisym- 
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Figure  4.  Evolution  history  of  the  {100},  {110},  and  {111}  pole  figures  at  the 
center  of  the  impacted  surface  during  the  Taylor-cylinder  impact 
experiment.  Pole  figures  are  shown  along  the  impact  direction;  this 
direction  was  previously  the  transverse  direction  in  the  rolled  Ta 
sheet  material 
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metric  manner,  the  material  behind  (away  from  the  impacted  surface)  retains  an  amount 
of  anisotropy  from  the  initial  conditions.  Thus,  the  initially  circular  surface  becomes  “oval- 
ized.”  Figure  5(b)  shows  a  view  of  the  cylinder  from  the  surface  away  from  the  impact 
site.  The  circular  top  of  the  cylinder  is  shown  as  a  reference  to  the  ovalized  footprint  at 
the  impact  site.  Maudlin  et  al.  [4]  reported  the  ratio  of  the  major  axis  to  the  minor  axis  of 
the  elliptical  foot  to  be  1.17.  Using  their  anisotropic  plasticity  model  to  reflect  the  initial 
material  anisotropy,  the  authors  predict  this  ratio  to  be  1.15,  by  accounting  for  the  evolution 
of  texture,  the  current  calculation  predicts  the  ratio  as  1.18. 


Figure  5.  Two  views  along  the  impact  direction  of  the  Taylor  cylinder; 

(a)  contours  of  effective  plastic  strain  in  the  ovalized  footprint,  and 

(b)  the  ovalized  footprint  relative  to  the  initially  circular  cross 
section  which  is  retained  away  from  the  impact  surface 

In  this  case,  it  is  the  initial  anisotropy  of  the  material  away  from  the  impacted  surface 
that  adds  asymmetric  confinement  to  deformation  in  the  foot.  Thus,  from  a  structural 
analysis  point  of  view,  the  texture  evolution,  although  extreme  on  the  impact  surface,  is  not 
a  major  contributor  to  the  ovaled  footprint. 

4.3  Explosive  Deep-Drawing 

In  this  section,  the  explosive  deep-drawing  of  a  cup-like  shape  from  an  initially  flat  disk 
is  discussed.  This  is  a  highly  simplified  version  of  many  other  explosive-forming  geometries, 
yet  relevant  in  that  we  can  examine  through-thickness  thinning,  and  thus  the  onset  of  both 
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diffuse  and  localized  failure  modes  as  a  function  of  the  initial  arrangement  of  Ta  single 
crystals  within  the  macroscale  aggregate. 

Initially,  a  thin  (relative  to  the  specimen  radius)  disk  of  pure  Ta  is  simply  set  on  top  of  a 
high-explosive  cylinder.  The  explosive  is  initiated  on  the  cylinder  face  opposite  the  Ta  disk, 
and  the  resulting  transmitted  shock  imparts  both  a  macroscopic  acceleration  to  the  disk  and 
deformation  into  a  cup-like  configuration.  To  enhance  the  deep-drawability  of  the  final  cup 
shape,  it  is  highly  desirable  for  the  workpiece  to  be  resistant  to  through-thickness  thinning. 
The  initial  geometry  is  represented  in  the  bottom  of  Figure  6  by  an  axis3nnmetric  finite 
element  mesh.  The  Ta  disk  is  represented  by  240  symmetric  crossed-triangular  elements 
with  12  elements  through  the  thickness.  The  initial  grain  distribution,  shown  in  Figure  7, 
is  assumed  to  be  random  and  homogeneously  distributed  throughout  the  disk.  At  50  /xs 
after  initiation  of  the  high  explosive,  the  original  disk  shape  is  deformed  into  a  cup  with  an 
upward  velocity  of  approximately  760  m/s.  Contours  of  effective  plastic  strain  (top  portion 
of  Figure  6)  show  that  Viigb  strain  regions  after  deformation  are  both  in  the  nose  and  tail 
regions  of  the  deformed  cup. 

After  deformation,  individual  grains  within  the  Ta  microstructure  begin  to  accumulate 
around  certain  preferred  orientations  (i.e.,  we  begin  to  develop  texture).  As  the  texture 
evolves,  it  also  becomes  graded  along  the  radial  direction.  The  center  row  of  Figure  7  shows 
the  evolved  texture  in  the  nose  region  where  a  strong  {111}  fiber  texture  has  developed  due 
to  the  predominant  uniaxial-compression  deformation  in  this  region.  This  is  consistent  with 
our  previous  result  and  also  with  those  of  Lee  et  al.  [9],  who  show  that  Ta  under  uniaxial 
compression  evolves  both  a  strong  {111}  fiber  texture  and  a  weak  {100}  fiber  texture.  In 
the  current  calculation,  it  is  the  small  number  of  grains  that  prevents  us  from  picking  up  the 
weak  {100}  component.  Note,  however,  that  5  grains  out  of  48  are  present  within  10®  of  the 
center  in  the  {100}  pole  figure.  Also  note  that  stress  calculations  take  place  at  the  element 
centroid,  off  the  axis  of  symmetry,  and  thus  deviate  from  ideal,  uniaxial  compression. 

The  bottom  row  of  Figure  7  is  taken  firom  the  high-strain  tail  region  depicted  in  Figure  6. 
In  these  pole  figures,  we  see  the  presence  of  a  strong  {111}  fiber  and  a  weaker  {100}  texture 
except  that  both  are  rotated  away  from  the  the  projectile  axis  due  to  large  shear  strains  that 
develop  in  this  region. 

The  previous  example  is  illustrative  of  the  kind  of  texture  evolution  that  can  be  expected 
in  Ta  during  explosive  loading.  Of  more  direct  importance,  however,  are  the  effects  of  initial 


17 


A=  3.00E-01 


Nose  region 


High-strain  tail  region 


time  =  50  |is 


Ta  disk 


High  explosive 


time  =  0 


Figure  6.  The  explosive  deep-drawing  process;  an  initial  finite  element  mesh 

(bottom)  depicting  a  Ta  disk  and  high-explosive  charge.  Initiation  of 
the  explosive  imparts  both  an  upward  velocity  to  the  disk  and 
deformation  into  a  cup-like  configuration.  Efiective  strain  contours  in 
the  explosively  deep-drawn  disk  at  50  fis  are  shown  at  the  top  of  the 
figure 
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Nose  region 
time  =  50  MS 


Figure  7.  Texture  evolution  during  explosive  deep-drawing.  Columns  show  the 
{100},  {110},  and  {111}  pole  figures.  The  top  row  characterizes  the 
initially  random  distribution  of  grains  that  is  homogeneous 
throughout  the  disk.  By  50  /zs,  there  is  a  gradient  in  texture  from  the 
nose  (center  row)  extending  radially  outward  in  the  deep-drawn  cup. 
Texture  in  the  high-strain  tail  region  is  shown  in  the  bottom  row 
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orientations,  and  thus  material  anisotropy,  in  mitigating  material  failure  during  this  kind  of 
forming.  In  a  series  of  articles,  Clark  et  al.  [31,  32,  33]  examine  the  effects  of  thermomechan¬ 
ical  processing  variables  on  the  production  of  Ta  sheet  products.  Based  on  observations  that 
bcc  steels  [34]  with  a  strong  {111}  fiber  orientation  have  enhanced  deformation  character¬ 
istics  during  mechanical  deep-drawing,  the  authors  examined  combinations  of  upset  forging 
and  cross  rolling,  as  well  as  various  heat  treatments,  with  the  aim  of  producing  a  strong 
{111}  fiber-oriented  Ta  sheet.  The  authors  examined  only  processing  variables  on  the  sheet 
product;  they  did  not  discuss  the  deep-drawability  of  these  products.  They  were,  however, 
able  to  show  that  different  dominant  peaks  in  the  orientation  distribution  function  (ODF) 
for  the  sheet  could  be  tailored  by  varying  the  processing  schedule. 

Following  the  work  of  these  authors,  it  is  desirable  to  understand  the  effects  of  various 
initial  textures  on  the  explosive  deep-drawing  characteristics  of  Ta  sheet.  Perfect  {100}, 
{110},  and  {111}  fiber  textures  were  produced  by  aligning  a  grain  with  its  {100},  {110},  or 
{111}  plane,  respectively,  normal  to  the  plane  of  the  sheet.  Random  rotations  within  that 
plane  produced  the  initial  distributions  of  48  grains  shown  in  the  first  column  of  Figure  8. 
These  initial  fiber  textures  were  incorporated  into  the  previous  amalysis.  It  is  felt  that  by 
understanding  how  each  individual  texture  component  manifests  itself  in  terms  of  macro¬ 
scopic  flow  and  resistance  to  thinning,  we  can  gain  a  better  understanding  of  the  optimal 
microstructure  for  this  application. 

The  second  column  in  Figure  8  shows  how  the  initially  ideal  fiber  textures  evolve  in 
the  apex  of  the  explosively  formed  cup.  Both  the  {100}  and  the  {111}  fibers  remain  intact 
throughout  the  deformation.  Since  the  predominant  deformation  in  this  region  is  uniaxial 
compression,  this  is  in  direct  agreement  with  Lee  et  al.  [9]  who  observe  that  both  {111}  and 
{100}  fiber  textures  are  stable  under  uniaxial  compression.  Lee  et  al.  [9]  also  observe  that 
{110}  fiber  orientations  evolve  into  {111}  fibers  under  uniaxial  compression;  this  can  be  seen 
in  Figure  8,  where,  at  50  /js,  the  {110}  fiber  in  the  nose  region  is  shown  to  be  in  the  process 
of  evolving  into  a  {111}  fiber.  Grains  originally  oriented  w33‘’  from  the  center  of  the  pole 
figure  are  now  moving  toward  the  center.  Further  uniaxial  deformation  would  complete  the 
transition  into  an  ideal  {111}  fiber.  In  the  tail  region  (last  column  of  Figure  8),  the  effects 
of  a  higher  shear-strain  environment  can  be  seen  on  the  initial  fiber  textures. 

Again,  it  is  not  really  these  textures  or  their  final  state  per  se  that  is  macroscopically 
important,  but  their  effects  on  the  drawing  process  itself.  Figure  9  shows  that  after  only  3 
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Figure  8.  {Ill}  pole  figures  showing  texture  evolution  due  to  three  different 
ideal  fiber  textures.  The  initial  homogeneous  textures  (first  column) 
correspond  to  ideal  {100},  {110},  and  {111}  fiber  textures.  The 
second  column  shows  final  deformation  textures  in  the  nose  region 
while  the  third  column  shows  final  texture  in  the  high-strain  tail 
region 
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fis  into  the  calculation,  the  through-thickness  stress  levels  in  the  {111}  and  {110}  oriented 
disks  are  significantly  higher  than  the  disk  with  an  initial  {100}  fiber  orientation.  These 
stress  levels  directly  refiect  relative  material  anisotropy  due  to  the  different  fiber  textures 
and,  consequently,  the  resistance  to  through-thickness  thinning  during  the  process. 

As-6.00E+0S 


initial  {100} 
fiber  texture 


initial  (110) 
fiber  texture 


initial  (111) 
fiber  texture 


Figure  9.  Contours  of  deviatoric  through-thickness  stress  (units  are  in  pascals) 
plotted  on  the  deformed  configuration  of  the  explosively  formed 
material  at  3  /xs 

Figure  10  shows  effective  strain  contours  on  the  deformed  shapes  once  formation  is  com¬ 
plete.  By  relative  comparison,  the  workpiece  with  the  initial  {100}  fiber  texture  underwent 
much  greater  through-thickness  deformations  than  either  of  the  {110}  or  the  {111}  fiber- 
oriented  workpieces.  Specifically,  large  through-thickness  thinning  is  apparent  throughout 
the  workpiece,  with  localizations  occurring  in  both  the  nose  and  tail  region.  The  {110}  and 
{111}  fibers  produced  approximately  the  same  gross  thinning  through  the  cross  section  with 
the  exception  that  the  {110}  fiber  produced  a  cup  showing  strong  tendency  for  localization 
within  the  nose  region,  while  the  {111}  fiber  gave  only  a  slightly  stronger  diffuse  necking 
region  in  the  tail.  It  is  also  interesting  to  compare  the  forming  result  from  the  {111}  fiber- 
oriented  material  to  that  of  the  random-oriented  material.  Figure  11  reproduces  these  two 
calculations  on  opposing  sides  of  the  symmetric  plane.  It  can  be  seen  that  the  randomly 
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oriented  material  undergoes  higher  through-thickness  strains;  in  particular,  the  nose  region 
is  approximately  30%  thinner  than  the  nose  region  in  the  {111}  fiber-oriented  material. 
Overall,  of  the  four  calculations  (three  fiber  textures  and  one  randomly  oriented  material) 
the  {111}  fiber  seems  to  be  the  microscale  arrangement  of  choice  for  explosive  deep-drawing 
of  bcc  materials. 


5.  Discussion 

The  current  article  establishes  a  framework  for  modeling  the  high-strain-rate,  nonhomo- 
geneous  evolution  of  anisotropy  in  polycrystalline  Ta  that  is  applicable  to  a  broad  class  of 
ductile  materials  that  deform  primarily  by  bulk  dislocation  motion.  Beginning  with  a  rate- 
and  temperature-dependent  model  for  dislocations  overcoming  both  intrinsic  lattice  resis¬ 
tance  and  barriers  due  to  dislocation  accumulation,  the  model  focused  on  shear-stress/shear- 
strain  behavior  for  slip  systems  in  Ta  single  crystals.  Incorporation  of  the  restricted-glide 
kinematics  and  an  appropriate  numerical  method  for  single  crystals  allowed  for  the  resolu¬ 
tion  of  individual  crystal  rotations,  and  thus  the  evolution  of  stress  anisotropy,  within  the 
polycrystal. 

Applications  of  the  model  were  designed  (1)  to  validate  the  model  by  way  of  connection 
to  other  work  (since  there  is  very  little  work  on  high-rate  texture  evolution  and  mechanical 
anisotropy  in  pure  Ta,  model  validation  is  not  a  straightforward  task);  (2)  to  show  that  the 
model  is  capable  of  resolving  nonhomogeneous  evolution  of  texture;  and,  most  importantly, 
(3)  to  illuminate  anisotropic  material  behavior  in  Ta  that  is  directly  relevant  to  modern 
engineering  applications. 

For  simple  compression,  it  was  shown  that  the  current  model  depicted  texture  evolution 
consistent  with  Lee  et  al.  [9].  The  cylinder  impact  example  was  chosen  to  show  that  the 
model,  while  reflecting  the  expected  texture  evolution  from  a  rolling-type  texture  to  that 
which  evolves  under  uniaxial  compression,  also  exhibits  macroscopic  phenomena  distinct  to 
anisotropic  materials  (i.e.,  the  ovalization  of  an  initially  circular  footprint).  The  extent  of 
ovalization  was  consistent  with  experimental  work  [4].  Finally,  a  detailed  examination  of 
an  explosive-forming  process  revealed  sensitivity  of  the  final  deformed  shape  to  the  initial 
microstructure.  In  particular,  it  was  shown  that  an  initial  {111}  fiber  texture  resulted  in 
a  more  stable  through-thickness  configuration  than  either  the  {100}  fiber  texture  or  the 


23 


As  3.00E-01 


B=  4.13E^r 


Figure  10.  Effective  strain  contours  on  the  deformed  configuration  for 

explosive  deep-drawn  material  at  50  /xs.  The  initial  textures  were 
distributed  homogeneously  in  the  disk  and  then  allowed  to  evolve 
with  deformation. 


FnmOE  UEVELS 


Initial  random-oriented  Initial  {111}  fiber-oriented 

material  material 

Figure  11.  A  comparison  of  effective  strain  contours  over  the  deformed  shapes 
from  the  initially  random-oriented  material  and  the  material  with 
an  initial  {111}  fiber  texture.  The  two  calculations  are  shown  on 
opposite  sides  of  the  plane  of  symmetry  for  the  calculation 

initially  random  microstructure.  The  relevance  here  is  twofold.  Clark  et  al.  [31,  32,  33] 
showed  that  under  common  stock  processing  conditions,  the  extent  of  {111} <  uvw  >  versus 
{100}<  uvw  >  orientations  is  a  process- dependent  material  property.  Therefore,  efforts  to 
bias  the  initial  material  processing  toward  creating  a  strong  and  homogeneous  {111}  fiber 
texture  are  not  in  vain.  Further,  since  the  initial  {111}  fiber  texture  resulted  in  better  flow 
properties  than  even  the  initial  random-oriented  material,  current  efforts  to  produce  stock 
Ta  via  powder  metallurgy  methods  may  result  in  a  suboptimal  product. 
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